***********************************************************
** Wrap-up paper: Estado e desigualdade
* Fluxos + previd�ncia + diferencial p�blico privado na POF
***********************************************************

***********************************************************
* Caminhos
	* Diret�rio de trabalho
		global dir "d:\bases_trabalhos\stateineq_2012"
	* Bancos de dados
		global dados "d:\bases_trabalhos\stateineq_2012\arq_para_analises.dta"
	* Arquivo auxiliar
		global aux "D:\Dropbox\trabalhos\Ipea\2012\15. Artigo - State and inequality\stata\aux_analises.do"
	* Vari�veis
		* Renda recebida
			global bruto "funcionalismo fp_simulado fp_dif fp_difpos fp_difneg seg_desemp previbruta previnss prevpubl ppmenteto ppmaiteto sim_ate sim_mais progsociais trabalho outrasrendas estadobruto rendabruta routras"
		* Tributos e contribui��es
			global tribs "estadotributos tributos previtribs cp_inss cp_rpps cp_rpps_ate cp_rpps_mais cp_rpps_sim_ate cp_rpps_sim_mais "
		* Renda l�quida
			global liquido "previliq inssliq rppsliq estadoliquido rendaliquida"
		* Para a decomposi��o: mais desagregado
			global maisdesagrega1 fp_simulado fp_difpos  fp_difneg seg_desemp previnss ppmenteto sim_ate sim_mais progsociais trabalho outrasrendas tributos cp_inss cp_rpps_ate cp_rpps_sim_ate cp_rpps_sim_mais
			global maisdesagrega2 fp_simulado fp_difpos  fp_difneg seg_desemp previnss ppmenteto ppmaiteto progsociais trabalho outrasrendas tributos cp_inss cp_rpps_ate cp_rpps_mais
		* Para a decomposi��o: m�dio desagregado
			global medio1desagrega fp_simulado fp_dif seg_desemp  previnss prevpubl progsociais routras tributos cp_inss cp_rpps
			global medio2desagrega fp_simulado fp_dif seg_desemp  previbruta progsociais routras tributos previtribs
			global medio3desagrega funcionalismo seg_desemp inssliq rppsliq progsociais routras tributos 
			global medio4desagrega funcionalismo seg_desemp previliq progsociais routras tributos 
		* Para a decomposi��o: pouco desagregado
			global poucodesagrega estadobruto estadotributos routras
		* Para a decomposi��o: nada desagregado
			global nadadesagrega estadoliquido routras
	* X% mais pobres e Y% mais ricos
		scalar xpob=50
		scalar yric=96
***********************************************************	


***********************************************************
* #0: Arquivo auxiliar
	qui do "$aux"
***********************************************************


***********************************************************
* #1: carrega dados (d'oh!)
	use "$dados", clear
***********************************************************


***********************************************************
* #2: M�dias gerais, UCs que recebem cada tipo e m�dias condicionais
	* m�dias gerais para todas as vari�veis
		qui medias $bruto $tribs $liquido [iw=pesopessoasuc]
		matrix define output2 = e(medias)
	* % que � afetado por cada vari�vel
		qui rec  $bruto $tribs $liquido [iw=pesopessoasuc]
		matrix define output2 = output2, e(rec)
	* m�dias condicionais (s� para quem recebe)
		qui medcond  $bruto $tribs $liquido [iw=pesopessoasuc]
		matrix define output2 = output2, e(medcond)
***********************************************************


***********************************************************
* #3: M�dias gerais, UCs que recebem cada tipo e m�dias condicionais
*		para os X% mais pobres e para os Y% mais ricos 

	* A) 50% mais pobres
		* m�dias gerais para todas as vari�veis
			qui medias $bruto $tribs $liquido [iw=pesopessoasuc] if cem<=xpob
			matrix define output3a = e(medias)
		* % que � afetado por cada vari�vel
			qui rec  $bruto $tribs $liquido [iw=pesopessoasuc] if cem<=xpob
			matrix define output3a = output3a, e(rec)
		* m�dias condicionais (s� para quem recebe)
			qui medcond  $bruto $tribs $liquido [iw=pesopessoasuc] if cem<=xpob
			matrix define output3a = output3a, e(medcond)
		
	* B) 5% mais ricos
		* m�dias gerais para todas as vari�veis
			qui medias $bruto $tribs $liquido [iw=pesopessoasuc] if cem>=yric
			matrix define output3b = e(medias)
		* % que � afetado por cada vari�vel
			qui rec  $bruto $tribs $liquido [iw=pesopessoasuc] if cem>=yric
			matrix define output3b = output3b, e(rec)
		* m�dias condicionais (s� para quem recebe)
			qui medcond  $bruto $tribs $liquido [iw=pesopessoasuc] if cem>=yric
			matrix define output3b = output3b, e(medcond)
***********************************************************


***********************************************************
* #4: Output dos resultados anteriores
	log using "$dir\output_medias", replace
	
	* Toda a popula��o
		matrix list output2, f(%10.2f)
	* X% mais pobres
		matrix list output3a, f(%10.2f)
	* Y% mais ricos
		matrix list output3b, f(%10.2f)	
	
	log close
***********************************************************


***********************************************************
* #5: decomposi��o do gini
	* LOG
		log using "$dir\output_ginideco", replace
	* Mais desagregada 1, 2
		ginideco $maisdesagrega1 [iw=pesopessoasuc]
		ginideco $maisdesagrega2 [iw=pesopessoasuc]		
	* M�dio desagregado 1, 2, 3 e 4
		ginideco $medio1desagrega [iw=pesopessoasuc]
		ginideco $medio2desagrega [iw=pesopessoasuc]
		ginideco $medio3desagrega [iw=pesopessoasuc]
		ginideco $medio4desagrega [iw=pesopessoasuc]
	* Pouco desagregado
		ginideco $poucodesagrega [iw=pesopessoasuc]
	* Nada desagregado
		ginideco $nadadesagrega [iw=pesopessoasuc]
	* Log close
		log close
***********************************************************		
